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Evaporatively cooled air ventilation for airconditioning (ECAVA) is an energy ef- 
ficient environmental friendly method as it uses only water as the working fluid. 
Though ECAVA is not a substitute of mechanical air conditioning, yet it renders 
quite a reasonable level of comfort for relieving the thermal strain of human beings 
and that too without using ozone depleting chemicals. 

A mathematical model developed by Verma was modified for calculations of the 
heat gaine .. by a building of a given geometry. Thereafter the ventilation air is 
calculated for a preassigned temperature rise of evaporatively cooled air. The 
salient feature of the present mathematical model is that it has provisions for 
making various changes in parameters through sets of input commands. The model 
was applied for a three storey production shop of a company at Kanpur involved in 
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manufacture of scooters and other components. The volume flow rates of air were 
calculated for all the three floors. For the known volume of air for the building 
the sizes of axial-flow fans were chosen from the existing catalogues of standard 
manufacturers. To reduce energy requirements and to have less vibration and 
noise the low speed axiaHlow fan was identified. The electrtic energy consumption 
was further reduced by introducing humidity controller. The same not only saves 
electric energy, but also reduces wear and tear of pump, heating of motor and 
noise drastically. Solenoid valves were also incorporated to help conserve electric 
energy. 

Since there are no laid down norms in India specifying the details of density for 
wood fibres used in evaporative coolers, the ASHRE guideline was followed. To 
keep the cooler pads always wet water distribution was ascertained by conducting 
experiments for the water flow rates to help the designer obtain practical ranges 
of discharge coefficients. 

The air distribution ducts were designed on the basis of equal friction pressure 
drop method. Care has been taken to balance the pressure in main branches 
as far as possible. To calculate various duct dimensions, a computer programme 
was developed. Wherever structural constructions came, the duct dimensions were 
changed maintaining the same diameter of the ducts. This approach was adopted 
to ensure equal friction drop for the system. 
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Chapter 1 
Introduction 


1.1 Description 


According to the Oxford Dictionary the art of ventilating is, “to expose to 
fresh air", and to, “cause air to circulate freely in an enclosed space”. The 
term general ventilation has acquired a special meaning in contradiction to 
‘local exhaust ventilation’, and refers to the commonly understood process in 
which the space is constantly flushed with fresh or humidified out door air. 

The object of air distribution in a ventilation and air conditioning system is to 
create the proper combination of temperature, humidity and air motion in the 
conditioned room for obtaining comfort conditions. The standard limit for an 
effective draft temperature’, has been established, for this purpose. The same 
consists of air temperature, air motion, relative humidity and there physiolog- 
ical effects on the human body[l]. Any variations in the accepted standards 
of one of these elements or lack of uniform conditions in the same pant of the 
space as room air temperature variations (horozontally, vertically or both) and 
excessive air motion (draft) will lead to discomfort. 

In India, Malhotra [2] has obtained the effective temperatures for hot & humid 
as well as hot & dry climates based on recording of 75 % votes by subjects 
about their sensation in different working environments. Table 1.1 illustrates 
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the comfort range found by him 


Table 1.1: Recommended Effective Temperatures [°C]. 


SI. No. 

Level of Comfort 

Effective Temperature 



Hot humid climate 

Hot dry climate 

1. 

Warm and unpleasant 

27.0 

26.0 - 28.3 

2. 

Comfortable and pleasant 
(upper level) 

24.0 - 25.0 

24.4 - 26.6 

3. 

Comfortable and pleasant 
(lower level) 

22.0 - 22.5 

21.1 - 24.3 


Table 1.1 shows that high effective temperature of 27 °C for subjects of hot- 
humid climate as compared to about 22 °C of existing practice. The energy re- 
quirement is reduced significantly due to the higher temperature. The thermal 
comfort requirements for hot-humid and hot-dry climates have been studied 
elaborately by Kimura and Tanabe[3] with a view to energy conservation in 
buildings. The results are presented in Table 1.2 having various combinations 
of human comfort parameters. The higher is the inside condition, the lower is 
the thermal potential across the building structure. Hence energy' requirement 
gets reduced considerably. 


Table 1.2: Parameters for Human Comfort [3]. 


SI. No. 

Experimental Condition 

Y elocity(m/s) 

1 . 

27 °C, 4> = 50 % 

0.5 

2. 

29 °C, <j> = 50 % 

1.2 

3. 

31 °C, <£ = 50 % 

1.6 


If high comfort condition is permitted, even evaporative cooling becomes fea- 
sible. In fact it was this very concept due to which the concept of hybrid air 
conditioning was envisaged. The hybrid system is a process based air condi- 
tioning system as the comfort condition in a hot-dry climate is achieved by 
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evaporative cooling rather than the wrongly used mechanical airconditioner. 
The mechanical airconditioner is not only the wrong system for maintenance 
of comfort condition in hot-dry climate, but also takes energy about 7 times 
that of an evaporative cooling system [4]. 

It is a well known fact that the evaporative airconditioning is not a substitute 
of mechanical airconditioning, but it renders reasonable degree of comfort and 
helps relieve thermal strain to a great extent. That is why these days, venti- 
lation with evaporative cooling is becoming popular to increase productivity 
in industrial organisations. 

For maintenance of desirable temperature, the circulation should be high enough 
to have less temperature rise. For this purpose total heat gain of the confined 
space or building is calculated. Thereafter the ventilation air is estimated on 
the basis of energy balance between the heat gain of the building and cooling 
to be produced by evaporativelv cooled air. 

The cooling heat gain calculations are based on the design outside maximum 
and minimum temperatures. However several researchers have represented 
outside daily hourly temperature by sinusoidal function. On the other hand 
the ASHRAE suggests an hourly temperature as the percentage of daily range 
for cooling heat gain calculations [5]. It gives the hourly multiplying factor 
in a tabular form for calculation by the designer. The ASHRAE approach 
is similar to using sinusoidal temperature variations based on maximum and 
minimum values. Because in both methods only maximum and minimum daily 
temperatures are needed and they dispense with the complete temperature- 
time history of the place for which heat gain is to be calculated. The hourly 
temperatures have been calculated using the method of ASHRAE and sinu- 
soidal temperature variations approach. In either method the predicted values 
differ by only a small magnitude about 2 - 3 %. The predictions for a few 
hours by the latter is compensated by the over prediction by the same for the 
rest of the period. The same is not going to make much changes in the final 
predictions. 
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1.2 Energy Conservation And Global Warm- 
ing 

The era of industrial revolution has arroused a need to be inquisitive as to 
why so many governments in the world have focused together to address the 
problem of global warming? The answers are not absolutely clear, but it is 
apparent that mankind’s activities are perturbing the natural variations of 
carbon, nitrogen and chlorine compounds in the atmosphere and many gov- 
ernments especially in the developed w T orld, have decided to take action in a 
big way. 

In late 1973 the Organisation of Petroleum Exporting Countries (OPEC) ex- 
tended the life of their income source by restricting oil production and allowing 
demand to force the price in excess of $ 2.50 per barrel for the major importing 
countries. This action led rapidly to revaluation of all forms of energy. The 
developed world focussed first on turning of lights and equipment when not 
required, drastically cutting back ventilation air quantities, and other similar 
conservation measures. 

It was not until a second, a much more substantial increase in price of oil 
in 1978, that the developed nations resorted to design and development of 
efficient goods and services. Evaporative airconditioning technologies even 
though they are still underutilised and often even unknown in many parts of 
the world, provide comfort cooling throughout many arid zones of the world, as 
well as relief cooling for commercial and industrial applications such as green 
houses, laundries, warehouses, factories and poultry houses[9]. 


1.3 Three Important “E’s”. 


. Thev stand for : 
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• Energy ( First E) -.Energy demand rising to sustain various industrial ac- 
tivities and other applications. 

• Environment Second E) increasing energy use renders environment de- 
gradation below allowed limit due to release of global warming and pol- 
luting gases, etc. 

• EconomyiThhh E) investment renders recovery to sustain existing need 
and future requirement. 

A few decades back the main aim was to meet the energy needs as far as pos- 
sible. Hence, the first E dominated. However, the increasing energy- use caused 
tremendous increase in release of particulate and polluting gases in addition to 
global warming gases. In India alone the carbon dioxide release from thermal 
power plants turns out to be over 800 million tonnes per year and in the world 
about 7 billion tonnes per year. Four main gases: nitrogen, oxygen, argon and 
carbon dioxide form about 99.99% of the atmosphere, with traces of rare gases 
like neon, methane, krypton, helium, xenon, hydrogen, carbon monoxide, radon 
and ozone. Of these gases COi , water vapour and methane are the green house 
gases (GHG). These gases allow the incoming shortwave radiation from the sun 
to the earth’s surface while the GHG’s absorb partly 7 the reradiated low en- 
ergy radiation, thereby 7 the environmental temperature near the earths surface 
remains at an average temperature of 288K as against 259K which it would 
have remained. Since the GHG’s surpassed the accepted limits of environmen- 
tal quality and the degraded environmental quality started affecting health of 
living beings including plants and crops. Hence, energy demand for sustaining 
various activities cannot be seen in isolation of pollution, globalwarming and 
ozone depletion. The ash released from the thermal power plants also cause 
pollution as well [6]. 

The pollution level has gone so high that the environmentalists are taking 
action against the reckless growth of thermal power plants. The seriousness of 
the environmental problems has also been accepted by the world bank. The 
same is reflected from the fact that the World Bank keeps the environmental 
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aspect mandatory for release of funds in support of development of new power 
plants. Thus, the second E stands for environmental pollution. 

Refrigeration systems cause two types of contribution to the global warming 
effect : 

• Direct Global Warming Potential (DGWP). 

• Indirect Global Warming Potential (IGWP). 

1. DGWP: Taking the reference base as the GWP of CO 2 over a 100 years 
period: 

(a) The GWP of Rl34a is 1200 times to that of CO 2 . 

(b) The GWP of R404A is 3500 times to that of CCT 

2. Total Equivalent Global Warming Impact (TEWI):The Total 
Equivalent Global Warming Impact includes both GWP, under condi- 
tions of refrigerant loss and recovery, expressed as ( GWPx total emis- 
sion). and the “additional” (so called indirect) global warming from CO 2 
produced throughout the system’s service life by generating the elec- 
tricity, or the direct combustion, needed to run the equipment. TEW I 
is useful for showing whether the direct effect or the indirect effect is 
dominant in a particular system [9] . 

TEW I is calculated from the equation: 

TEWI = {GWP) x (M) + (a) x (B) (1.1) 


where, 

GWP = GWP of refrigerant relative to CO 2 and GWP {CO 2 ) = 1 
M = Total mass of refrigerant released (kg). 

a = Amount of CO 2 released in generating electricity ( kg CO 2 per KW) 
B = energy consumption of the system in its life time. 
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TEW I is much more reliable indicator of global warming for a particular 
system than is GW P. 

CFC refrigerants being responsible for ozone depletion/ozone hole. The 
Evaporative Cooling Airconditioning (EGA) system requires only wa- 
ter as the working medium. Hence, the system does not release global 
warming gas like CO 2 nor polluting gases. Even this system takes only 
one-fifth to one-eighth of the electric power as compared to the vapour 
compression system . Hence, the systems TEWI is very less for ECA. 
When water-vapour compression system is used, in this case also green 
house gas emissions are absent. But, water-vapour itself is a global warm- 
ing gas. 

With the increasing competetiveness in business, to compete with the 
other quicker methods of cooling, the investment in the evaporative cool- 
ing system renders much recovery to sustain future requirement. Hence, 
evaporatively cooled air for airconditioning (ECAA) has acquired great 
signifiance as the investment for power for ECAA is much smaller than 
that for the vapour-compression system rendering saving in revenue for 
further upkeep, upgradation and instillation of new power plants to meet 
increasing demands in future. This aspect is very important. Evapo- 
rative cooling systems are more inexpensive to purchase and easier to 
operate as compared to the vapour-compression airconditioners. 

1.4 Present Work 

A computer programme is modified for calculation of cooling loads for a 
manufacturing shop in the Industrial Estate, Kanpur. The unsteady na- 
ture of the outside temperature variation and the incident solar radiation 
is included by considering a decrement factor and time lag factor. By 
this approach the calculations are simplified significantly. The details are 
available in Sec.2.3. Based on the heat gain for the building the volume 
of ventilation air required per floor of the three storeyed manufacturing 
shop was worked out, keping the space limitations and user’s needs in 
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mind. Depending on the volume of air flow per floor a computer pro- 
gramme for duct sizing was developed keeping the distribution of air to 
all parts of the building in view. 

Finally, evaporative cooling systems were designed and their layout was 
suggested for instillation on top of the building due to space restrictions. 
Experiments were conducted for evaluation of discharge coefficient for the 
holes in the trough used for water distribution to the evaporative coolers. 
To keep energy conserved and without affecting evaporative cooler per- 
formance electronic timers have been recommended for all three evapora- 
tive cooling systems. The air handling capacities for evaporative cooling 
systems for ground, first and top floors turned out to be, 80,000m 3 /h, 
40,000m 3 /h and 14,400m 3 /h, respectively. 



Chapter 2 

Cooling Load Calculation 


In this chapter.the expressions used for evaluation of the various individ- 
ual heat gains are first presented followed by a summation of all loads to 
give the total heat gained by the structure and finally the expression for 
the volume of ventilation air are presented. 

2.1 Different Heat Sources : 

The cooling load for the evaporative air conditioning systems for summer, 
in general, comprises: 

• Sensible Heat Load 

• Latent Heat Load 

The following constitute the sensible heat load [10] : 

(a) Solar Radiation : It is again divided into two categories: 

i. Direct Solar Radiation : It is the heat load which results from 
the transmission of solar radiation through glass windows and 
ventilators. 

ii. Indirect Solar Radiation : A part of the solar radiation falling 
over the walls and roofs gets absorbed and then transmitted to 
the room by conduction and convection. 
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(b) Heat Transfer by Temperature Difference: The heat trans- 
fer through the exterior walls, roofs, floors, doors etc., takes place 
due to temperature difference between the surroundings and interior 
space. 

(c) Heat Released from Occupants 

(d) The electrical energy dissipation from various equipments and ap- 
pliances add to the heat energy. 

(e) Miscellaneous heat sources. 

The sources which contribute to the latent-heat load are listed below: 

(a) The latent heat load from the occupants; 

(b) Moisture passing directly into the space through permeable walls 
where water- vapour pressure is high; 

(c) The latent-heat load from cooking foods and from stored materials. 
The cooling load arises from the following sources: 

(a) External Sources : 

i. Heat transmission through barriers such as walls, doors, ceiling, 
floor, etc. being caused by the temperature difference existing 
on the two sides of the barrier. 

ii. The solar heat , absorbed by walls , roofs and windows exposed 
to radiation from the sun and transferred to the inside space. 

iii. . Heat and moisture introduced to the conditioned space through 

infiltration and ventilation air. 

(b) Internal Sources : 

i. Occupancy load, both sensible and latent 

ii. Lighting load 

iii. Power equipment load, e.g., fans, water pumps, etc. 

iv. Appliances load, both sensible and latent 
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2.2 Need for Hourly Cooling Load Calcu- 
lations 

It is a well known fact that the ambient air temperature is not con- 
stant during day, and accordingly the cooling load varies over a period 
of 24 hours. This is mainly caused by the variation in the solar intensity 
falling on the earth surface. The cooling heat gain calculation is further 
complicated by the fact that a wall has thermal capacity, due to which 
a certain amount of heat passing through it is stored and is transmit- 
ted to the inside at some time latter. Therefore, calculation based on 
instantaneous heat transmission through structure without considering 
the thermal capacity of the wall, does not reflect the actual heat gained 
by the confined space. Moreover, the traditional methods of evaluating 
the cooling load involves various assumptions like the load from each 
component is constant in a day and maximum value of the cooling load 
is equal to the individual maximum, which, in general, is not correct, as 
the occurrence of maximum heat gain at particular time is not the occur- 
rences of the other components. It is proper to adopt an hourly cooling 
load calculation to facilitate an economic selection of an evaporative air 
conditioning unit [5]. 

2.3 Structural Cooling Load : 

2.3.1 Building Survey 

The building survey incorporates the flowing particulars: 

(a) Location of the building ,i.e., longitude and latitude 

(b) Orientation 

(c) Dimension of the building structure, such as length, breadth , height 
and thickness of each layer of the building material 

(d) Composition of building material and their physical properties. 
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2.3.2 Hourly Outside Temperature Variation 

Figures 2.1 to 2.4 show the hourly variation in outside temperature for 
Kanpur. The meteorological data reveals that the minimum temperature 
occurs just one or two hours before the sunshine while the maximum 
temperature occurs three to four hours after the solar noon [15]. Since, 
the temperature time history is available for only a few places, an hourly 
variation in temperature is predicted based on maximum and minimum 
temperatures. 

The ASHRAE Handbook of Fundamentals [5] has outlined a procedure 
to predict the hourly temperature. In the present work, the same is 
followed. 


t 0 — T max — Daily Range x P ercentageF actor / 100 (2.1) 

where, f 0 =temperature at any time in °C 
DailyRange — Umax TAm 

The percentage factor is taken for that hour for which outside temper- 
ature is to be calculated. These factors are tabulated in [5] and are 
reproduced here in Appendix ‘A’. 

The comparison of the temperature variation by theoretical with actual 
variation is shown in figures 2.1 to 2.4. The data have been taken from 
the average temperature of three years, i.e., 1982, 1983 and 1985 [15], 
and are reproduced here in Appendix : B’. 

2.3.3 Solar Heat Gain 

Solar radiation forms a significant part of the cooling load for a building. 
The total radiation (I t ), reaching on terrestrial surface is the sum of the 
direct solar radiation (Id), the diffused sky radiation (Id) and the solar 
radiation reflected from the surrounding surface (I r ). The intensity of 
the direct component is the product of the direct normal radiation (Idn) 
and the cosine of the incidence angle ( 8 ) between the incoming solar rays 
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Figure 2.1: Hourly variation of actual and theoretical temperature in April. 
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Figure 2.3: Hourly variation of actual and theoretical temperature in June. 

and a line normal to the surface as shown in figure 2.5. Thus, 

It = ID+Id + Ir ( 2 - 2 ) 

In the present analysis the reflected radiation (J r ) is neglected. 

Id is given as: 

Id = -Tow x cos9 
where, Idn is expressed as : 

^ DN exp(B/ sinfi) 

Id consists two parts: 

(a) diffused solar radiation from clear sky 

Ids = C.Idn-Fss 


(2.3) 

(2.4) 


(2.5) 
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time,h 

Figure 2.4: Hourly variation of actual and theoretical temperature in July. 






2.3 Structural Cooling Load : 


16 


where. F ss is the angle factor between the surface and the sky. i.e., 
the fraction of short wave radiation emitted by the sky that reaches 
the tilted surface (dimensionless). F ss is 0.5 for vertical surfaces 
and 1.0 for horizontal surfaces. For other surfaces, 

F S s = (1 + cosS)/2 (2.6) 

where, E is the tilt angle from the horizontal plane(Fig 2.5). 

(b) ground reflected diffused solar radiation 

Idg ~ J-tH -Qg-Fsg ( 2-0 

where, g 9 is the reflectance of foreground and F sg is the angle factor 
between the surface and the ground (dimensionless), and F sg = 
1 —F ss . I t H represents the total solar radiation falling on the ground, 
and is given asl*# = Ids{C + sin 3). 

The value of A, B and C are taken from [5] and the same are repro- 
duced in Appendix ‘C’ for the twenty first day of the each month. 


2.3.4 Determination of Angles of Incidence 


The sun’s position in sky is most conveniently expressed in terms of solar 
altitude 8, above the horizontal and solar azimuth angle <p, measured 
from south. These angles 3 and p are expressed in the terms of the 
latitude of the place l, solar declination d, and hour angle h. 


sinP = cosl.cosd.cosh + sinl.sind 


and, 

, sinfi.sinl — sind 
COS9= cotp.cosr 
The hour angle (in degree) is calculated as : 


h = (12 — LST) x 15 


( 2 . 8 ) 

(2.9) 


( 2 . 10 ) 


If a surface is tilted by an angle E to the horizontal plane, then the 
incidence angle 0 is given by 


cosd = cosp.cosj.sinT, + sinfi.cosll 


( 2 . 11 ) 
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When the surface is horizontal. E = 0°. and : 

cosd H = sinB ( 2 . 12 ) ?£ 

For a vertical surface, E = 90°, and : 

cos9y = cosp.cosj (2.13) 

The surface solar azimuth angle 7 is given as : 

7 = cj) — ib , before solar noon (2.14) 

7 = <fi + ip , after solar noon (2.15) 


where, ip is surface azimuth angle which represents the orientation of 
the surface. The ip values for different orientation are ips( 0°), b’£(90°), 
v.v(180°) and ^(270°). 

2.3.5 Problem Formulation 

2. 3. 5.1 Heat Transfer through “Walls and Roof 

Heat Transmission through the wall or roof of building structures is not 
steady and is therefore, difficult to evaluate. The two principal factors 
causing this are: 

(a) The variation in the outside air temperature over a period of 24 
hours. 

(b) The variation in the solar radiation intensity that is incident upon 
the surface over a period of 24 hours. 

The phenomenon is further complicated by the fact that a wall or roof 
has a thermal capacity due to which a certain amount of heat passing 
through it is stored and is transmitted to outside or inside at some later 
time. 

The problem requires a solution of the governing equation for unsteady- 
state one-dimensional heat transfer, viz., 

dt d 2 t 

3r ~ Q X dx 2 


(2.16) 
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where, t is the temperature at any section of the wall or roof at a distance 
x from the surface at a time r, and a is the thermal diffusivity given by 


a = 


k_ 

p.c 


(2.17) 


where, k is the thermal conductivity and p.c is the heat capacity of 
the wall or roof, in which p and c are the density and specific heat, 
respectively. 

The heat transfer equation is to be solved with the boundary conditions 
of periodic variation of outside air temperature and solar radiation. 


2.3.6 Sol-air Temperature 

The calculation of structural load as a result of heat gain through exterior 
roof and walls involves the concept of sol-air temperature [6]. A heat 
balance at a sunlit surface gives the heat flux to the surface as: 


q = a s I t + h 0 (T o - T so ) - eAR (2.18) 

e = emissivity of the surface, 1.0 for black body 
A R = 63.0 W/m 2 , for horizontal surface and 0, for vertical surface [5] 
a s = absorptivity of the surface Equation 2.17 can also be written as 

q = ho(T sol - T so ) (2.19) 


then, from equations 2.17 and 2.18 : 

Tsoi = T 0 + a s It/h 0 - e ARJh a 

Thus, the combined effect of the solar radiation and outside air temper- 
ature has been incorporated into a single effective temperature known as 
sol-air temperature. 

2.3.7 Heat Transfer Coefficients for Walls and Roofs 

The outside film coefficients are given by [8,16]: 
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For exterior walls 


h 0 = 7.373 + 5.066 x V 0 , ( W/m 2 .°C ) 

(2.20) 

For roof 



h Q = 7.9953 + 6.364 x V 0 , ( W/m 2 .°C ) 

(2.21) 

where, V 0 is the outside air velocity in m/s. 

The inside film coefficients are given by [8]: 

For interior walls 



hi = 1.77 x (At) 0 ’ 25 + (7.373 + 5.066 x l/) 

,(ir/m 2 .°C) 

(2.22) 

For roof 



hi = 1.31 x (At) 0 - 25 + (7.9953 + 6.364 x 1/) 

, (W/m 2 .°C) 

(2.23) 


where, 

At — | t s { t{\ , 

t si = inside surface temperature, 

U = inside temperature, and 
'll = inside air velocity in m/s. 

2.3.8 Heat Transfer Through Walls And Roofs 

If the thermal capacity of the wall is ignored, then the instantenous rate 
of heat transfer through the wall at any time is given by 

Q = UA(T e - Ti ) (2.24) 

However, most building materials have a finite thermal capacity C which 
is expressed as 

C = me = pcV — pc(AAX) (2.25) 

where, 

m = Mass of wall, [kg] 
p = Density of wall material, [m 3 /kg] 
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.4 = cross sectional area of wall material. [m 2 ] and 
AX = wall thickness [m]. 

It has been seen that there is a two-fold effect of the thermal capacity 
on heat transfer 


• a time lag between the heat transfer at the outside surface and at 
the inside surface 

• decrement in heat transfer due to the absorption of heat by the wall 
and subsequent transfer of a part of this heat back to the outside 
air when the temperature of outside air is lower [10]. 


when these facts are taken into account, the actual heat transfer through 
the structure at any instant ‘t’ is given by the following relation[2.26] 


Qt = B- 4 A(?e - Ti) + AjUjXjiTo - T e )} 

j= 1 


(2.26) 


where, 

T e = mean sol air temperature °C, 

A = decrement factor, 

T 0 = temperature at time 0 before the time ; t’, 
o = time-lag (h). j = j th surface. 

If thickness Ax is less than 0.30m then 

, 0.8 x Ax 


(2.27) 


or A = 0.1 

U , A, and T 0 refer to the overall heat transfer coefficient, cross sectional 
area and temperature at time t — (f>, respectively. The summation symbol 
X) stands for the four walls and the roof. 

When the wall is very thin, the time lag factor A j — > 1 and we get 


Q t = j:A j U j X j (T 0 ~T i ) 


i=i 


(2.28) 
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On the other hand, in case of very thick wall A j — > 0 and the heat transfer 
is obtained from. 


Q i = Y.A j U 1 (T,-T i ) 


7—1 


(2.29) 


The overall heat transfer coefficient U is given by 


U 


1 3 Ax ■ 

_ x , y" —Zi 
h ^ Z-' k- 

jz=l 


1 

hi 


(2.30) 


T 0 = ambient temperature(°C) Tj = inside temperature(°C) Axj and kj 
= thickness and thermal conductivity of the ‘jth’ layer of the structural 
material hi = convective heat transfer co-efficient for inside surface. kW/mVC 
U = overall heat transfer co-efficient, kW/m 2 .°C Figures 2.6 and 2.7 give 
the structural detail of walls and roof, respectively. 


CEMENT PLASTER 



A x2 = 0.115 m 


kl = 0.00072 

k2 = 0.00130 
lc3 = 0.00072 


x3 = 0.012 m 


Figure 2.6: Structural detail of wall. 
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MORTAT 



kl = 0.00072 


MORTAR 

k2 = 0.00173 
k3 = 0.00072 

Figure 2.7: Structural detail of ceiling. 

2.3.9 Heat Transfer through Glass 

Glass construction forms a significant part of modern building structures. 
The heat transfer through glass [10] comprises: 

• all the transmitted radiation 

• a part of the absorbed radiation that enters the conditioned space, 
and 

• the heat transmitted due to temperature difference between the 
outside and inside temperatures. 

The direct radiation enters the space only if the glass receives the direct 
rays of the sun. The diffuse radiation enters the space even when the 
glass is not facing the sun. Figure 2.8 represents all sorts of heat transfer 
into conditioned space through glass. The heat transfer to the space is 
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give by : 


Q glass A sun . i p.Ifc) ‘ ~r hi.A w i n .(tgi ti) 


(2.31) 


where, 

tgi = temperature of inner surface of glass 
A sun = glass area directly exposed to the sun 
A u ~i n = total glass area of window 

The subscripts D and d denote the terms for direct and diffuse radiations, 
respectively. 



Writing energy balance for glass sheet itself, we have 

A sun .ccD.lD d - A w i n .a.d.Id = A w i n \hi .{t gi t,) ~r h 0 .[tg O to)] (2.32) 
where, i ao = temperature of outer surface glass 
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As the thickness of glass sheet is very small, inner and outer surface 
temperatures of glass can be assumed to be equal to (f s ), without creat- 
ing appreciable error in cooling load but simultaneously simplifying the 
analysis tremendously. Now, modified equation can be expressed as : 

, A st m-OiD-^D ■^w 'in • Oid-Id + Awin-h'i-ti T ^win-h 0 .t 0 ^ 

l g _ — ■» - [2.66) 

■s^win \'h » iL o) 


Putting the value of t 9 in place of t 9 i in equation (2.45) ,we get the cooling 
load due to glass as : 


r) a ~ T i A _ r t Asuti-QD’Id + Awin'&d-Id , TT A fj. j. ^ 

glass — ^sun-^D-ijD * ' d^d * “ cU *A. win .\t 0 t{ ) 


I 1 - + fr) 

where, U is the overall coefficient of heat transfer given by 


(2.34) 


h 0 h{ 


(2.33) 


If the thermal resistance of glass is also considered, U is given by : 


_1 _ _1_ Ax 1 
U ~ h~ 0 + ~k~ + 'h i 


(2.36) 


where, Ax is the thickness of glass and k g is its thermal conductivity. 

Here, it has to be kept in mind that the transmissivity and absorptivity 
are function of the angle of incidence. They vary for direct radiation 
while are almost constant for diffused part. These values are given in 
Appendix ‘D’[13]. 


2.3.10 Shading of Surfaces 

The most effective way to reduce the solar load on fenestration is to 
intercept direct radiation from the sun before it reaches the glass. To 
serve this purpose, most glass areas are provided with reveals, overhangs 
and fins in the form of vertical and horizontal projections from the walls. 

The ability of these projections to intercept the direct component of 
solar radiation depend on their geometry, surface-solar angle (7) and the 
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profile or shadow-line angle (fi) (Fig 2.9). The profile angle is defined 
as the angular difference between a horizontal plane and a plane tilted 
about a horizontal axis in the plane of the fenestration until it includes 
the sun. The profile angle can be calculated by: 


tanQ = 


tanB 
cos 7 


(2.37) 


The shadow height. Sh, due to horizontal projection. P, on a window 



Figure 2.9: Shading on glass due to horizontal and vertical projections. 


or w r all for any time of day or year is related to the profile angle 9. by: 

S H = P x cotQ (2.38) 

The shadow width S w , due to a vertical projection. P r . (it need not not 
be equal to the horizontal projection) on a window or wall for any given 
time of day and year is related to the vertical surface -solar azimuth angle 
9' by: 


S w = P v x cotj 


(2.39) 
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Thus, the sunlit area of the window is: 

= (W - S U .).(H - S H ) (2.40) 

where, W and H are the width and height of window. 

2.3.11 Heat Transfer through Door 

The cooling heat gain due to doors are given by 

Qdoor = U-A-door-iLo tj) (2.41) 

where, U is the overall heat transfer coefficient through door [5] and A door 
is area of the door. 


2.4 Ventilation Load 


Ventilation air is mandatory to ensure fresh air supply in the residential 
buildings. The amount of ventilation needed has been debated for over 
a century, and the different types of rationale developed have led to 
radically different ventilation standards. The current rationale [5] for 
the minimum outside air requirement is 2.5 L/s per person based on 
CO 2 concentration. The commonly accepted value is 2.364 L/s. 

Ventilation, also, causes both types of load namely sensible heat gain 
(Qvs) and latent heat gain (Q v i). The load due to ventilation are given 
by following expressions : 

Q vs — n x 1.232 xLxAT , (TV) (2.42) 


and, 


Q vl = nx 3012 x L x Acu , (TV) 


(2.43) 


where, AT is difference between outside and inside temperatures, °C. Aw 
is difference between inside and outside air specific humidities, kgw./kg d.a.. 
L is ventilation air in litre per second per person and n is number of peo- 
ple in the conditioned space. 
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2.5 Occupancy Load 

The occupants in a conditioned space give out heat at a metabolic rate 
that depends on their rate of working. The relative proportion of the 
sensible and latent heats given out. however, depends on the ambient 
dry bulb temperature. The lower the dry bulb temperature, the greater 
the heat given out as sensible heat. The instantaneous sensible load is 
the product of the sensible heat loss from the people and the cooling load 
factor (CLF). This CLF [5] is function of the time people spend in the 
conditioned space and the time elapsed since first entering. The sensible 
cooling load (Q os ) and latent load ( Qoi ) can be expressed as : 

Q os = n x Sensible heat loss x CLF , (W) (2-44) 

and, 

Qoi — n x Latent heat loss , (W) (2-45) 

where, n is number of people in the conditioned space. 

2.6 Lighting Load 

Electric lights generate a sensible heat equal to the electric power con- 
sumed. Most of the energy is generated as heat and the rest as light 
which also becomes heat after multiple reflections. An accurate esti- 
mate of the cooling load imposed by the lighting is not straight forward. 
Because the rate of heat gain to the air caused by lights can be quite 
different from the power supplied to the light points or lighting fixtures. 

The time lag effect should be taken into account in calculating the cooling 
load, since the actual load is lower than the instantaneous heat gain, and 
peak load may be significantly affected. 

The lighting load (Qis) is expressed as : 

Qis = Light Wattage x Special Allowance Factor x CLF (2-46) 
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where, cooling load factor (CLF) [5] is a function of time of use , type of 
arrangement, room furnishing, etc. Special allowance factor is introduced 
for fluorescent fixtures and fixtures requiring more energy than their 
rated Wattage. This is due to the fact that the choke takes about 20 - 
25% of the rated power of the tube, used in mercury vapour tube lights. 


2.7 Power Equipment Load 

Power equipment load ( Q ps ) is calculated from 

i _ kW Rating x Load Factor x 1000 x CLF 
ps %Motor Efficiency / 100 

where, load factor is merely the fraction rated load delivered under the 
conditions of the cooling load estimate. 

Cooling load factor (CLF) [5] is the function of time (hours after equip- 
ments are on) and total operational time. 


2.8 Appliances Load 

Most appliances contribute sensible and latent heats. The latent heat 
produced depends on the functions the appliances perform. Electric mo- 
tors contribute sensible heat to the conditioned space. A part of power 
input is directly converted into heat due to the inefficiency of the motor 
and is dissipated through the frame of the motor. The motor efficiency 
of 80% is considered in the case study presented. In estimating cooling 
load, heat gain from heat producing appliances , e.g. kitchen appliances, 
computers, etc., must be taken into account. Some appliances produce 
only sensible heat load while others generate sensible (Q as ) as well as 
latent load (<Q q i). 

Qas = Sensible Heat Rate x CLF (2.48) 


and, 


Qai = Latent Heat Rate 


(2.49) 
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Here, it can be noticed that sensible load is present even after putting 
off the appliance while latent heat load is present only for that hours in 
which appliance is on. 


2.9 Total Cooling Load 


The solar load (Q ss ) is the summation of 

Qss = Q structural ~ Qglass d" Qdocrr (2.50) 

Total cooling load is sum of all the sensible loads and latent loads. 

Q sensible ~ Qss d~ Qvs d" Qos ~ Qls d~ Qps d~ Qas (2.51) 


and, 


Therefore, 


Qlatent — Qvl 1 Qal ' Qol 
Qtotal ~ Q sensible ~ Qlatent 


The design cooling heat load is S.F. x Q to tai where, 
S.F. = safety factor =about 1.1. 


(2.52) 


(2.53) 


2.10 Ventilation Air Calculations 


The mass flow rates of ventilation air is calculated from the total cooling 
load Q for the conditioned space. Thus the evaporatively cooled ventila- 
tion air requirement for a temperature difference of AT is calculated by 
[16] in the following expressions: 


Q = m a x Cp x AT , (klV) 

Q 


m„ — 


( kg/s ) 


Cp x AT 5 
V = m a v ,(m 3 /h) 


(2.54) 

(2.55) 

(2.56) 



2.10 Ventilation Air Calculations 
where, 

?77 n =mass flow rate of ventilation air requirement, (kg/s) 
V'=Yolume flow rate of ventilation air required , (m 3 /hi) 
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Chapter 3 

Evaporative Airconditioning 
And Duct Design 


Air cooling by water evaporation occurs in nature near waterfalls, flowing 
air streams over lakes and oceans, under summer showers and even upon 
wetted skin. The evaporative process simply removes sensible heat (i.e., 
cooling by decreasing the surface temperature) and replaces it with latent 
heat (i.e., increasing the moisture content of air). It is interesting to note 
that the sea acts as an evaporative cooler due to which the temperature 
in costal areas of India hardly goes beyond 30 - 35°C as compared to 40 - 
45°C in the hot-dry regions. This is supported from the detailed results 
presented in [18]. 

3.1 Theory 

Evaporation is described as an adiabatic process, during which the to- 
tal thermal energy contents of the system remains constant. Figure 3.1 
shows the depression in temperature resulting from latent heat gain of 
air during evaporation of water. As water evaporates, the sensible heat 
content of the system falls, while the latent heat content increases by an 
equal amount. In other words, dry-bulb temperature of air falls, but its 
moisture content rises. The limit of temperature reduction is up to the 
wet-bulb temperature of the air at the beginning of the process. Surface 
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DRY AIR 


MOIST AIR 
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WATER 
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KEAT 
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WATER 


YOU FEEL 35*C 
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LATENT HEAT ENERGY 


YOU FEEL 21°C 


Figure 3.1: Schematic representation of depression in dry-bulb temperature 
during evaporative cooling [7]. 


evaporation under ideal conditions can cool up to the wet-bulb temper- 
ature. The process stops when the relative humidity of air approaches 
100 %. A simple measure of the potential for evaporative cooling at any 
given air condition is the wet-bulb depression, defined as the difference 
between the dry-bulb and wet-bulb temperatures. It provides an upper 
limit of the achievable temperature drop. 

Water is a universal coolant and the working fluid that can meet air- 
conditioning needs in both residential, commercial and industrial ap- 
plications. Two principle methods of evaporative airconditioning are 
commonly used: 

• Direct Cooling Process 

• Indirect Cooling Process 

Direct cooling in which water evaporates directly into the air stream, 
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thus reducing the air dry-bulb temperature, while humidifying the air. 

Indirect cooling, where primary air is cooled sensibly with a heat ex- 
changer, while the secondary air carries away the heat energy from the 
primary air as generated vapour. Direct and Indirect processes can be 
combined (indirect/direct). Compared to vapour compression systems, 
increased air flow rates are used for direct evaporative comfort cooling 
to compensate for higher supply air temperatures. 

3.1.1 Direct Evaporative Airconditioning 

In direct evaporative cooling air is drawn through wetted pads or a spray 
chamber and its sensible heat energy causes evaporation of water, the 
dry-bulb temperature of air gets reduced. The air temperature is de- 
creased by 60% to 90% of the wet-bulb depression (ambient dry- bulb 
temperature less wet-bulb temperature.) Note that there is no sensi- 
ble cooling and that this is essentially' an isenthalpic process for direct 
evaporative airconditioning (also termed evaporative, “cooling' ). 

In arid regions direct coolers provide comfort cooling, while high humid 
areas can use direct cooling for specialised applications. Direct EAC 
consumes significantly' less energy than vapour-compression refrigeration. 
The only power consuming components of an evaporative cooler are fans 
and small water pumps. Energy savings of evaporative coolers vary with 
humidity levels and temperature. Direct systems in low-humiditv zones 
typically realize an energy r saving of 60% to 80% over refrigerated sys- 
tems. In-direct evaporative cooling equipment, water is supplied through 
a float valve through a small reservoir from where it flows through fibrous 
pads. A fan draws large volume of outdoor air through the pads, where 
it is cooled by evaporation, and then supplied to the building. This cool 
and more humid air absorbs sensible heat from the buiding. An efficient 
wetted pad can reduce the air temperature by as much as 95% of the 
wet-bulb depression, while an inefficient and poorly designed pad may 
only reduce this by 50%, or less. 
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Figure 3.2: Schematic representation of direct evaporative cooling system. 

v 

In the actual evaporative cooling process there is very little change in the 
energy state of air. Direct ECA is simple and cheap due to the reduced 
energy requirement for fan and make up water, but it is ineffective if am- 
bient wet-bulb temperature is close to dry-bulb temperature, the cooling 
effect is not sufficient enough for indoor comfort cooling applications, but 
it finds some applications (e.g., green houses and industries spot cooling). 
Direct evaporative coolers should not recirculate indoor air as cooling ef- 
fectiveness decreases. 

The saturation effectiveness of a direct evaporative airconditioner best 
describes the performance of an ECA unit. Saturation effectiveness is 
described as the difference between the entering and exit dry-bulb tem- 
peratures over the -wet-bulb depression and can be defined as [7]: 

T] s& — T sl( ppiy)/H' Bdepression (3-1) 


where, 
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77 se =saturation efficiency 

^ Pdepression = (Pdb ~ Pub) outside 

r d6 =outdoor dry-bulb temperature 
T wb = outdoor wet-bulb temperature 

3.1.2 Indirect Evaporative Cooling 

In indirect evaporative cooling water does not evaporate in the air sup- 
plied into the conditioned space. It attempts to make use of the evap- 
orative cooling process without increasing the amount of moisture in 
the supplied air. Indirect evaporative cooling equipments use a heat ex- 
changer to avoid the direct contact between water to be evaporated and 
the supply air. A direct evaporative process cools air that flows across 
one side of the heat exchanger, removing heat, and is then exhausted to 
the atmosphere. The air to be supplied to the building flows across the 
other side of the heat exchanger and is cooled without receiving any 
moisture (Fig. 3.3) The processes are: 

Process 1 - 2 : Humidification of secondary air. 

Process 2 - 3 : Heat from primary air transfered to secondary air. 

Process 1 - 4 : Cooling of primary air for sopply to the conditioned room. 

The second direct stage is then added to further cool the air. In this 
system the outside air is pre-cooled in an indirect stage and then further 
cooled in a subsequent direct stage. The first stage cools the air without 
adding moisture and in second stage air passes directly into the wrnter 
spray system. This generally yields the final air temperature leaving 
the evaporative cooler about 3.5°C lower than what could be achieved 
with a direct ECA only. This expands the ECA considerably to areas 
with slightfy higher wet-bulb temperatures. Generally, effectiveness of an 
indirect stage evaporative cooling system is found to be 65% is reached 
which allows ambient wet-bulb temperature of upto 25°C to provide low 
enough room temperature for real comfort cooling. 

The performance of the indirect evaporative cooling is measured by PF 



3.1 Theory 


36 



DRY - BULB TEMPERATURE oC 


HUMIDITY 

RATIO 

Kg MOISTYRE 
PER Kg DRY AIR 


Figure 3.3: Schematic representation of different processes taking place during 
evaporative cooling 


which is the ratio of the reduction of the dry-bulb temperature of the 
dry side air stream to the initial difference between the dry' side drv-bulb 
and wet side wet-bulb temperatures [7]. 

PF — (Tdb — T 'supply) /kb B depression (3-2) 


where, 

PF =performance factorncy 

11 B depression — ( Tdb ~ Twb)outside 

Tdi=outdoor dry-bulb temperature 
2 n stip pi y =supply air dry-bulb temperature 

An indirect evaporative cooling process is illustrated in fig. 3.5 In this 
method the return air at a temperature lower than the ambient temper- 
ature is used to cool the fresh air in a surface-to-surface heat exchanger. 
The system as well as the processes are shown in figures 3.4 and 3.5 
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Figure 3.4: Schematic representation of plate type Indirect/Direct evaporative 
cooling airconditioner 


respectively. 

The processes are: 

Process 3 - 4 : Humidification of primary air from the room. 

Process 4 - 2 : Heat gained from the room by the humidified air 

Process 1 - 3 : Sensible cooling of primary air by the secondary return 
air from the room. 

Process 2 - 5 : Sensible heat gain by the room return air in the surface 
to surface heat exchanger HE1. 

By this arrangement the temperature of the humidified air is further 
lowered by a few degrees. Hence, it provides better comfort. If it is 
desired to get an even lower temperature, the return air can be passed 
through a h umidifi er. By this approach the dry-bulb temperature of the 
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Figure 3.5: Schematic representation of indirect/direct evaporative cooling 
process 


return air is further lowered. As such it renders further reduction in the 
dry-bulb temperature of the primary air passing the surface-to- surface 
heat exchanger. Thus after humidification of this cooled air renders 
temperature which may be very close to the mechanical airconditioning 
state. The processes for the latter modified system is represented by 1 - 
3 - 4 - 2 - 6 - 7 in fig.3.6. 


3.2 Air Duct Design 

The essential economics of an air transmission system is achieved by a 
proper balance between the initial or final cost for the given flow rate of 
air. The first cost is detirmined by the cost of the duct system which 
depends on the duct size. The operating cost is detirmined by the fan 
power consumption which depends on the friction pressure drop in the 
evaporative cooling and duct system. The pressure drop can be reduced 
by increasing the size of the ducts but it will increase the first cost. 
Hence, there should be a proper balance between duct size and operating 
cost. A few general rules are stated below which may form as a guideline 
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Figure 3.6: Schematic representation of an indirect/direct evaporative cooling 
process 


for the design of ducts. 

• Air should be conveyed as directly as possible to economize power, 
material and space. 

• Sudden changes in direction should be avoided. When bends are 
essential, turning vanes should be used to minimize the pressure 
losses. 

• Air velocities in ducts should be within permissible limits to keep 
the noise level within allowable limits. 

• Diverging sections chould be made gradual. The angle of divergence 
should not exceed 20°[12,13]. 

• Rectangular ducts should be made as nearly square as possible. 
However, keeping the dimensions as constant as possible also helps 
reduce the use of material and construstion cost. This will ensure 
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minimum duct surface, and hence cost for same air carrving capac- 
ity. An aspect ratio of 4:1 or less should be maintained. 

• Ducts should be made of smooth material such as galvanised iron 
(GI) or aluminium sheet. Whenever other material is used al- 
lowance should be made for roughness of material. 

3.2.1 Duct Design Procedure 

There are three common methods for sizing of ducts. They are: 

(a) Equal Friction Method : In the equal friction method, the frictional 
pressure drop per unit length of the duct is maintained constant 
throughout the duct system. The procedure is to select a suitable 
velocity in the main duct based on practical requirements. Thus by 
knowing the air flow rate and velocity in the main duct, the size and 
friction loss are detirmined from the frangers charts or by using the 
expression of frictional pressure drop. The frictional pressure drop 
is found from: 

~ = 0.02061 ’ 1 ' 845 /ci 4 ' 92 (3.3) 

A L 

where, 

V'=volume flow rate of air, (m 3 /s) 
d=equivalent duct diameter, (m) 

Ap=pressure drop equivalent to water head, (mm) 

This expression gives the relation among pressure drop, volume-flow 
rate and equivalent diameter of duct. 

If the volume-flow rate is given, and assuming the allowed velocity 
of the duct, the diameter of the circular duct can be found from: 

d = x V/tv x v, (m) (3.4) 


where, 

u=velocity of air in the duct, (m/s) 

For known V, and d, the pressure drop per unit length is calculated. 
Thereafter, for any volume flow rate and known pressure drop d is 
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calculated as: 


d = ((Ap/AL)/[0.0206y 1 - 845 ]) 1 / 4 - 92 , (m) (3.5) 

This method of sizing ducts automatically reduces the air velocity in 
the direction of flow. The method is generally recomended because 
of simplicity. If an equal friction design has a mixture of short and 
long runs of ducts, the shortest duct will need considerable amount 
of dampering. This is a draw back of the equal friction design. 

(b) Velocity Reduction Method: In this method the main duct is de- 
signed in the same manner as in the equal friction method. There 
after arbitary reductions are made in the air velocities as we go 
down the duct run. The equivalent diameters are found as in the 
first case from the friction chart. Though the method allows safe 
velocities, it is not normally addopted unless the person using it has 
considerable practical experience and knowledge to design within 
reasonable accuracy. 

(c) Static Regain Method: The principle of the static regain method 
is to maintain a constant static pressure before each terminal and 
each branch. This is achieved by sizing the duct in such a manner 
that after each branch or outlet, the static pressure gain due to 
the reduction in velocity exactly balances the pressure drop in the 
succeeding duct design section [12, 13]. 


3.3 Rectangular Equivalent of Circular Ducts 

Air ducts are usually sized first for round sections. Then if rectangular 
ducts are required, ducts are designed to provide the same flow rates 
and to have the same rate of pressure drop as for round ducts. The 
rectangular equivalent of a round duct is given by[10,12,13] as: 


i? '= 1 - 3021 (iw r ’ (m) 


(3.6) 
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where, 

D £ =equi valent diameter of the circular duct,(m) 
a and 6=sides of the rectangular duct,(m) 

The equation may be used to detirmine one dimension of a rectangular 
duct if the other is assumed equivalent to a circular duct whose diameter 
is known. Although round ducts require the least metal to carry a given 
quantity- of air, rectangular ducts are used often because of the following 
reasons: 

• Space considerations as they fit easily in the building construction 
and occupy less building space without being conspicuous. 

• Ease of fabrication. 

Square ducts are closest to round ducts. They require less material than 
rectangular ducts. The material required for the rectangular ducts in- 
crease with the aspect ratio, viz a/b. 



Chapter 4 

Duct Design Procedure and 
Experimental Setup 

4.1 Duct Design and Layout 

The steps for calculation of the ventilation air are given in chapter 3 
sec. 3.6.for different volumes of ventilation air for each floor of the triple 
storey building, and the following procedure was adopted for the ducts: 

• Knowing the total volume of air required per floor, firstly air dis- 
tribution for the different parts of the building was planned on the 
basis of the users requirements. 

• The ducts layout was planned keeping the following as guidelines: 

— The space restrictions within the building. 

--Equal distribution of air for comfort of all occupants. 

— The users requirements. 

- Balanced system. 

• The friction charts of air were used to evaluate the equivalent di- 
ameter of the ducts keeping a constant pressure drop per length of 
duct and the known volume flow rate of the air as planned earlier. 

• A computer program was used to convert the equivalent diameter to 
a rectangular section using the expression given in chapter 3. The 


43 



4.2 Choics of Axial-Flow Fan 


44 


conversion charts were also used for this purpose. Care was taken 
to provide a nearly square section of the rectangular duct to incor- 
porate maximum saving of construction material but. at the same 
time due to space restrictions within the structure a compromise 
betweenw the two had to be resorted to. 

• All the bends were provided with curved vanes within the duct to 
reduce pressure losses within the ducts. 

• The air velocity in the main ducts were taken to be 9 m/s, being 
the allowable velocity for tolerable sound level for industrial appli- 
cations. Though the duct layout was found to be better and eco- 
nomical when the layout was as per the drawings attached. These 
figures were modified in order to satisfy the users due to the ma- 
chinery layout and lighting fixtures. These constraints cropped up 
after the design and drawing were completed. The modified air dis- 
tribution system ditributes air in every two bays and satisfies the 
users need. 


4.2 Choice of Axial-Flow Fan 

The catelogues of standard manufactures were used to choose the axial- 
flow fans. The following were incorporated while making the choice 

• Choice of low rpm fans to reduce noise levels. 

• Small diameter fans to facilate the space restrictions. 

• Delivery of the volume flow rate required per floor. 

• Low electric energy consumption. 

• Delivery of the required pressure head per floor. However in the 
present case the fan diameter is large for the required air delivery. 
The fan is a little costlier and the same is compensated by the 
reduced running bill. 
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4.3 Design of Evaporative Cooler System 

Based on various guidelines of ASHRAE Book of Equipment[20] for the 
range of parameters were designed for evaporative cooling systems. No 
strict norms are available, and the system design was made keeping the 
latest developments in mind. Generally straight evaporative pads are 
most commonly used having the water distribution at the top. The 
water distribution is done through the holes in the pipes. This type of 
system is followed for evaporative cooling systems installed in India. The 
experiments to arrive at few results which were needed were performed 
as explained subsequently. A wood wool pad used as an eliminator is 
of 30mm thickness and 0.9 kg/m 2 of the eliminator surface area. This 
is to prevent any water carry over to the supply duct. This becomes 
essential as in the existing evaporative cooling supply duct the carry 
over of water droplets to the supply duct and finally its leakage through 
the duct was a serious problem. This point is taken care of by introducing 
the eliminator. The eliminator thickness is reduced in order to reduce 
the pressure drop in the air supply. 

4.3.1 Experimental Setup 

A pump was used to circulate water from the collection tray through a 
calibrated rotameter to a trough of dimensions 1000 x 150 x 50. The 
trough was used to collect the water from the rotameter which was 
through a control valve. The water in the trough was then made to 
flow back to the collection tray through a number of holes of equal diam- 
eter drilled in the base of the trough.There were two rows of holes drilled 
in the trough 30mm apart. The distance between the holes in each row 
was kept at 60mm. The experiment was conducted for three set of hole 
diameter 2mm, 2.5mm, and 3.175inm respectively.The photographs of 
the experimental setup are attached below 




Figure 4.1: Experimental Setup 
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4.3.2 Experiment Details 

• Water was made to flow through the rotameter at different flow 
rates controlled by the adjustable valve. 

• For the different flow rates when the head of water in the trough 
became constant, the head and flow rates were measured and noted. 
This procedure was repeated for 2mm, 2.5mm, and 3.175mm diam- 
eter holes. 

In this maimer a set of readings were noted for the three different diam- 
eter hole sizes and a gcnral idea of how the Cd and head vary with the 
discharge. 

4.3.3 Density of Pads 

Since there are no laid down parameters for the density of wood wool to 
be kept per pad. A rough estimate given by ASHRAE is 1.5kg/m 2 to 
2kg/m 2 of surface area for a 50mm thick wood wool pad. 

In order to establishe the most appropriate pad density for the cooling 
system rccomended a frame of 30cm x 30cm x 5cm was prepared and 
varying densities of wood wool fibres were packed into the frame to choose 
the most appropriate density by judgement. The photographs of the 
same are attached for reference 




Chapter 5 

Results and Discussion 


Calculations for the present case study carried out for the manufacturing 
shop at Kanpur enable an experience as to how far there is difference 
between the general approach and the pratical situation. This chapter 
deals with results obtained by calculation of cooling loads on the basis of 
the theoretical and actual temperature variation of outside air. For this 
a generalised computer programme[ ] is modified in order to in corporate 
formulation for detirmination of ventilation air supply for different floors 
fo the triple storey building. 

5.1 Input Data 

the input data for the present problem are : building geometry di- 
mensions, orientation, location, machinery details, light fixtures, number 
of occupants etc. Hence, the total cooling load was found for different 
days and months. 

For calculation of the ventilation air, maximum of maxima found for 
different months was finally identified for which the ventilation air was 
detirmined. 
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Graphical Representation of Results 

5.2 Graphical Representation of Results 


5.2.1 Theoretical and Actual Temperature Varia- 
tion 

The actual hourly temperature variation as well as the suggested hourly 
variation in temperature based on maximum and minimum temperature 
variation of the day have been shown in figures 2.1 to 2.4. The theoretical 
temperature generally matches very well with the actual values except 
deviating here and there with 2 - 3 K. The deviation between them are 
such that over prediction at one time gets nullified by under prediction 
at other time. Moreover the present approach simplifies calculation pro- 
cedure significantly. 

5.2.2 Variation of Wet-Bulb Temperature 

Figure 5.1 shows the variation in wet-bulb temperature for the month of 
April, May, June and July for Kanpur city. There is not much variation 
in WBT in a particular day. The maximum WBT occurs at about 3 
P.M. except, in July, in which the maxima occurs at about 10 A.M. This 
variation may be different for different places. 

5.2.3 Variation of Relative Humidity 

Figure 5.2 shows variation in relative humidity with time for the month 
of April, May, June and July for Kanpur city. The relative humidity for 
these months is minimum at about 3 P.M. when dry-bulb temperature 
is maximum. 

5.2.4 Cooling Load Variations 

The total cooling load on the structure varies during the day as per 
magnitudes of the maximum and minimum outside temperatures. The 
variation of the total cooling load is represented by graph for each floor 
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Figure 5.1: Hourly variation of wet-bulb temperature. 

of the production shop as shown in figure 5.3. The maximum cooling 
load on the building occurs around 5 P.M., i.e.. about 5 hours after the 
peak solar load is felt outside. This is due to the heat capacity of the 
building and the various openings in the building. Figure 5.3 show the 
effect of the decrement factor and the time lag on the cooling load. It 
means that a building’s material and thickness of wall play a significant 
role in the solar heat gain by the building. It is also insteristing to note 
that though the instantaneous solar heat gain is very high, the wall ther- 
mal capacity helps reduce the peak cooling load, i.e., the lower capacity 
refrigeration equipment may be used to take care of cooling of the build- 
ing. Hence, the initial investment gets significantly reduced. But the 
s am e may not render significant reduction in the running cost, because 
the heat retained by the wall during the solar heat gain is reduced dur- 
ing the night, causing increased cooling compared to the existing lower 
thermal potential during the evening or night. 
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time,h 

Figure 5.2: Hourly variation of relative humidity. 

5.2.5 Experimental Results 

The graphs are plotted to show how the discharge coefficient varies for 
different hole diameters. A graph indicating head vs discharge is shown 
to ascertain suitable minimum head of water required to be maintained 
in the trough to keep the wood wool pads always wet. 

Experiments were conducted for three sizes of holes in the trough namely 
9 mm A 2.5mm^ and 3mm<£. The first size hole displayed a high head in 
the trough. At the same time the holes were getting blocked with sedi- 
ments in the water. On the other hand the head in the trough becomes 
very low for the 3.175mm<p holes. Thus, there is no uniform water dis- 
charge over the lenth of the trough. Hence, finally 33 numbers of 2.5mm<p 
hole size per meter was selected and recomended for instillation in the 
evaporative cooling system. The results are represented graphically in 
figure 5.4., and are also shown in table 5.4. 
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Figure 5.3: Variation of cooling load for May 

5.3 Tabulated Results 

The cooling loads for different floors of the triple storey building are 
given in Table 5.2 representing the maximum values. The electric load 
for each floor was taken as 40% of the connected load. The cooling load 
is found to be maximum in the month of May for which the air flow rates 
of the evaporatively cooled air required per floor of the building under 
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study were evaluated. 
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Table 5.1: Maximum Cooling Loads for Ground Floor and Total Load 


Month 

Structload 

Ventload 

Glassload 


Occuload 


April 

26.32 

3.06 

0.29 

94.96 

6.9 

131.53 

May 

38.08 

3.03 

0.59 

94.96 

6.9 

naEfiii 

June 

36.68 

0.00 

0.66 


6.9 

IBifclill 

July 

20.65 

1.74 

HUill 

94.96 

6.9 


Aug 

25.37 

7.17 




mi 


27.41 

6.20 

0.29 

94.96 

6.9 

135.77 

Oct 

12.64 

2.92 

0.04 

94.96 

6.9 

117.46 


Table 5.2: Maximum Cooling loads for First Floor and Total Load 


Month 

Structload 

Ventload 

Glassload 

Elecload 

Occuload 

Total 

April 

16.56 

4.08 

4.15 

18.19 

9.2 

52.90 

May 

26.05 

4.04 

4.36 

18.91 

9.2 

62.57 

June 

21.59 

0.42 

4.19 

18.91 

9.2 

57.31 

July 

12.29 

2.62 

3.78 

18.91 

9.2 

46.80 

mm 


mm& 

3.91 


9.2 

57.73 

Sept 

17.64 

— v\ i Mr 

3.46 

18.91 

9.2 

n 1 

Oct 

5.46 

3.89 

1.39 

18.91 

9.2 

1 


Table 5.3: Maximum Cooling loads for Second Floor and Total Load 



Structload 

Ventload 



Occuload 

Total 

H33PV 


1.79 



4.02 

98.33 

ipm 


1.77 

4.36 


4.02 

107.78 




2.91 

mmt 

4.02 

104.24 

i mm 


■na 

2.20 

mmi 

4.02 

90.97 

mm 


■n 

3.91 


4.02 

96.79 

Sept. 


3.62 

3.46 

1.46 

4.02 

95.46 

Oct 


1.70 

1.39 

1.46 

4.02 

73.79 


5.3.1 Experimantal Results 

The experimental results are tabulated for holes of different diameters. 
The measured flow rates for different heads of water in the trough are 
presented in Table 5.4 to 5.C. The discharge coefficient calculated from 
equations shows an increasing trend. The expression for Cd for different 
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diameters are given in equations. They can help calculate the discharge 
coefficient for the detirmination of the water required for circulation for 
the evaporative cooler pads. Hence, the total flow rate of water required 
for a given size of evaporative cooler pad can be estimated in this way. 

The expression for the coefficient of discharge is given by: 


• For 2mm <ji> hole size: 






C d = 4.3433 x 
For 2.5mm <}> hole size: 

C d = 2.7802 x 
For 3.175mm </> hole size: 

C d = 1.7237 x 


Q_ 

\[H 


Q 

VH 


Q 

JH 


where, 

Q — the discharge in U.S.gallons/min. 
H = the head in the trough in mm. 

C d = the coefficient of discharge. 


(5.1) 


(5.2) 


(5.3) 
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Table 5.4: Experimental results for 2mm 0 hole Size 


SI No 

Discharge(US Gal/min) 

Head (mm) 

Coeff.of Disch.(Cd) 

1. 

0.4 

8.5 

0.614 

2. 

0.5 

12.5 

0.614 

3. 

0.6 

18 

0.614 

4. 

0.7 

24 

0.620 

5. 

0.8 

30 

0.634 

6. 

0.9 

36 

0.651 

7. 

1.0 

41.5 

0.674 

8. 

1.1 

52 

0.662 

9. 

1.2 

57 

0.690 

10. 

1.4 

78.5 

0.686 

11. 

1.6 

105 

0.678 

12. 

1.8 

121 

0.717 

13. 

2.0 

130 

0.762 

14. 

2.1 

135 

0.78 


Table 5.5: Experimental results for 2.5mm 0 hole Size 


SI No 

Discharge(US Gal/min) 

Head (mm) 

Coeff.of Disch.(Cd) 

1 . 

0.4 

3.5 

0.594 

2. 

0.5 

6.0 

0.567 

3. 

0.6 

8.0 

0.589 

4. 

0.7 

10 

0.615 

5. 

0.8 

12.5 

0.629 

6. 

0.9 

' 15 

0.646 

7. 

1.0 

18.5 

0.646 

8. 

1.1 

20.5 

0.675 

9. 

1.2 

25.5 

0.660 

10. 

1.4 

33.5 

0.672 

11. 

1.6 

43.5 

0.674 

12. 

1.8 

54.5 

0.677 

13." 

2.0 

65.5 

0.087 

14. 

2.1 

74.5 

0.676 


5.4 Density of Wood Wool Pads 

The ASHR.AE Handbook of Equipment [7] gives that the density of 
wood fibers used for an evaporative cooling system should vary between 
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Table 5.6: Experimental results for 3.175mm 0 hole Size 


SI No 

Discharge(US Gal/min) 

Head (mm) 

Coeff.of Disch.(Cd) 1 

1 . 

1.2 

10 

0.654 1 

2. 

1.4 

13.5 

0.657 

3. 

1.6 

18 

0.650 

4. 

1.8 

21.5 

0.669 

5. 

1.9 

23 

0.682 

6. 

2.0 

25.5 

0.682 

7. 

2.1 

27.5 

0.690 

8. 

2.2 

30.5 

0.686 

9. 

2.3 

33.5 

0.684 

10. 

2.4 

37 

0.680 

11. 

2.5 

39.5 

0.685 


1.5kg/m 2 to 2.0kg/m 2 for a thickness of 5cm. of the pad surface area. 
The pad thickness should vary between 30mm and 55mm. In the present 
case 50min pad thickness has been recomended due to the large size of 
pad. An experimental result of reference [ 21 ] shows an optimum density 
of 1.97 kg/m 2 and 2.5kg/m 2 of the pad surface area. These values are 
even higher than those recomended by ASHRAE. Hence, it would render 
an increased pressure drop and hence reduce the air handling capacity of 
the fan. In the present case, however, the wood wool density of 1.5kg/m 2 
of the pad surface area based on the practical considerations and ex- 
periment is found to be appropriate for the evaporative cooling system 
recomended for the building. The photographs of the same are attached 
below. 


5.5 Axial-Flow Fan 

On reffering to the catelogues of standard manufactures the following 
axial-flow fans were recomended for instillation in the evaporative cooling 
systems of each floor. The head for fans in millimeters of water are 
recomended on the basis of the total pressure drop required for the air 
supply, evaporative cooling devices, ducts, etc. The details of the same 
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Figure 5.5: Photographs of recomended Wood Wool Pads 
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arc given below: 

• Ground Floor: 

(a) Discharge: 80,000 m 3 /h 

(b) Head: 30 mm of water. 

• First Floor: 

(a) Discharge: 40,000 m 3 /h 

(b) Head: 25 mm of water. 

• Second Floor: 

(a) Discharge: 14,400 m 3 /h 

(b) Head: 20 mm of water. 


5.6 Duct Design and Layout 

As explained in chapter4 the duct layout and size were calculated with 
the help of a computer program and the equal friction charts for air 
and GI sheets. The layout and sizes are as indicated in the drawings 
attached. 


5.7 Water Distribution 

Knowing all the parameters of the trough from the experiments con- 
ducted and refering to the catelogues of standard manufactures the fol- 
lowing was recomended 

• Ground Floor: 

(a) Rating: 0.75 kW 

(b) Discharge: 5000 L/h 

(c) Head : 6m of water 

• First Floor: 

(a) Rating: 0.37 kW 
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(b) Discharge: 2400 L/h 

(c) Head: 6m of water 

• Second Floor: 

(a) Rating: 0.18 kW 

(b) Discharge: 1200 L/h 

(c) Head: 6m of water 

5.8 Design of Wood Wool Pads 

• Ground Floor: The pads have two straight sections with curved 
pads at the back as given in Fig.5.6. 

• First Floor: Two straight pads each 2m long having a curved length 
of 3.3m length, each 2m high and 50mm thick having the mass of 
woodwool as 1.5kg/m 2 of the pad surface area. The details of the 
same are shown in Fig 5.7. 

• Second Floor: The pads have two straight pads 0.5m long having a 
semicircular pad with radius 0.75m, each being 1.5m high as shown 
in Fig 5.8. 

The evaporative cooling systems for all the three floors have been placed 
on the roof of the triple storey manufacturing shop due to the space re- 
strictions imposed within the factory. A detailed plan view and elevation 
of the systems are given in the drawings attached. It is worth mention- 
ing that the factory is located in an area where ammonia vapours from 
a fertilizer factory in the vicinity of the building under study causes sig- 
nificant nuisance. The evaporatively cooled air supply does away with 
this problem and provides a better quality of air for the users inside the 

workshop. 
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Door 600 X 1600 



Plan view 

Figure 5.6: Evaporative cooling system for ground floor 

5.9 Design of Water Distribution System 

A layout of the water distribution system proposed are given in the 
layout of the plan mew depicting the cooling systems placed on the roof 
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Elevotion 



Plan View 


Figure 5.7: Evaporative cooling system for first floor 

of the manufacturing shop. The water supply is divided along 3 paths 
for each evaporative cooler. Two water supply pipes for the straight 
pads and one for the curved pads. This arrangement helps operate the 
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Elevation 



Plan view 

Figure 5.8: Evaporative cooling system for second floor 
Figure 5.9: Plan view and elevation of the evaporative cooling systems 
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e.vaporat,ive system as usual with the water distribution to all pads in the 
summer season. But, in the rainy season when the humidity of ambient 
air is very high Kanpur (about 98%), the water distribution to the pads 
cum be restricted by stopping the water supply to any one of the pads. 

To get rid of continuous operation of the pump, an electronic timer is 
introduced. This arrangement renders the following advantages : 

• Wear and tear of pump reduces drastically. 

• Electric energy consumption reduces to about | t/l of the normal 
electric energy consumption. 

• Sound of the operating pump is reduced. 

• The motor runs always cool and hence the burning of the coil wind- 
ing in the motor is dispensed with. 



Chapter 6 


Conclusions and 
Suggestions 

6.1 Conclusions 

The following are the salient features of the present work : 

(a) A computer programme has been developed to predict the total 
cooling load on the basis of the hourly temperature variation 
based on T max and T mi „ for a given day. The vcrsality from 
the present approach is supported from the fact that only T max 
and T m i n values are needed in the computation and the same are 
available in several data books [5,16,17] for many places. But 
the hourly yariation is available only at selected meteorological 

centers. 

(b) The cooling load for the summer months is found to be max- 
imum in the month of May. This is in tune with the recom- 
mended month for practice engineers by the standard t,exts.[10,14,l 5] 
The ventilation air is eatimated on the basis of energy balance 
between the heat gained by the building and the cooling to be 
produced by the evaporatively cooled air. 

(d) The evaporative cooling systes are provided with eliminators 
between the pad and the axial-flow fans to avoid water droplets 
from entering the ducts and the conditioned space through the 
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ducts. 

(c) All bends, in the ducts have been provided with vanes to reduce 
picssure losses and dampers to control the air supply at desired 
rates. To reduce pressure losses a gradual expansion and con- 
traction is recomended for the duct whenever there is a change 
in section 

(f) The layout of the duct work was so designed that all occupants 
get the required volume of fresh evaporatively cooled air at the 
same time a simple network of ducts are proposed to help in its 
construction and instillation. 

(g) Solenoid valves have been provided in the water distribution 
system to help cnserve the wastage of water. The use of elec- 
tronic timers for the pumps used for make up water (manufac- 
tured by Mr Alok Prasad, at 1. 1. T. Kanpur) have been proposed 
which also help in conserving energy. 

(h) The straight pads and the curved pads are connected separately 
to water supply pipes so that in the summer season the evap- 
orative system will operate with the water distribution to all 
pads but in the rainy season when the humidity of ambient air 
is high, the water can be supplied to only one or two pads only. 
In this way the humidity can be controlled to some extent. 

(i) To avoid the instillation of an oversized pump for ground floor 
systems a stand by pump has also been recomended, which will 
be supplying water to the curved pads only. 

(j) The manufacturing shop is located in the industrial area in 
Kanpur and it was observed that there was a predominant odor 
of ammonia gas in the environment. The ammonia will be 
dcssolved in the water flowing down the evaporative cooler pads 
and hence, the system will supply better quality air to the work 
place. 
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6.2 Scope for Future Work 

(a) This computer programme to calculate the cooling load is ap- 
plicable to the buildings located in the Northern Hemisphere. 
This program can be suitably modified for the buildings located 
in the Southern Hemisphere by changing the subroutine which 
calculate the solar radiation intensity. 

(b) The problem can be extended to one used for all the year round 
air conditioning by the use of hybrid airconditioning system. 

(c) Economic analysis can be done for the evaporative cooling sys- 
tem to find initial investment and running cost, though it can 
be concluded that it will less than the cost of installing and 
running a central airconditioning system for the same building. 

(d) As discussed in chapter 1 more publicity should be given to the 
evaporatively cooled system because 

i. There is substantial energy and cost savings. 

ii. No chloroflurocarbon usage. 

iii. Reduced C0 2 , he. no ozone depletion, and power plant emis- 
sions, 

iv. Improved indoor air quality, 

v. It is easily intergrated in built up systems, 

vi. Potential for significant local fabrication and employment, 

(e) A more generalised programme needs to be developed which 

may design all ducts after taking input data. 
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Appendix A 

Percentage of the Daily 
Range 


Table A.l: Percentage of the Daily Range 


time,h 

% 

time,h 

% 

time,h 

% 


m 

1 


7 

93 

13 

11 

19 

34 

2 


8 

84 

14 

3 

20 

47 

B 


9 

71 

15 

0 

21 

58 

■Is 


10 

56 

16 

3 

22 

68 

5 

100 

11 

39 

17 

10 

23 

76 

G 

98 

12 

23 

18 

21 

24 

82 
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Appendix B 

Temperature and Relative 
Humidity of Outside Air 


Table B.l: Hourly Outside Temperature (°C) for Kanpur 
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Table B.2: Hourly Outside Relative Humidity (%) for Kanpur 






Appendix C 

Constants for Solar 
Radiation Calculation 


Table C.l: Constants for Solar Radiation Calculation 


Month 

Equation 

A 

B 

C 


of Time 


(Dimensionless 


min. 

W/m 2 

Ratios) 

Jan 

-11.2 

1230 

0.142 

0.058 

Feb 

-13.9 

1214 

0.144 

0.060 

Mar 

-7.5 

1185 

0.156 

0.071 

Apr 

1.1 

1135 

0.180 

0.097 

May 

3.3 

1103 

0.196 

0.121 

Jun 

-1.4 

1088 

0.205 

0.134 

July 

-G.2 

1085 

0.207 

0.136 

Aug 

-2.4 

1107 

0.201 

0.122 

Sep 

7.5 

1151 

0.177 

0.092 

Oct 

15.4 

1192 

0.160 

0.073 

Nov 

13.8 

1220 

0.149 

0.063 

Dec 

1.6 

1233 

0.142 

0.057 


,4 = Apparent Solar irradiation at zero air mass , (Wfm ) 
B = Atmospheric extinction coefficient 
C — Diffuse radiation factor 
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Appendix D 


Properties of Ordinary 
Glass 


Table D.l: Properties of Ordinary Glass for Direct Radiation 


Angle of Incidence 


Transmissivity, r D 

0.87 0.87 0.87 0.865 0.86 0.84 0.79 0.67 0.42 0.00 

Absoptivitv, np 

0.05 0.05 0.05 0.055 0.06 0.06 0.06 0.06 0.06 0.00 


For diffuse radiation, regardless of the angle of incidence 
r,t = 0.79 
a d = 0.06 
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Appendix E 


Input Data 

Table E.l: Input Data For Ground, First and Second Floors. 


| ShNo. ' 

Building Parameter 

SSlijSiil! 


Second Floor 


1. 

lat, long (°) 

nflii 

-do- 

-do- 

2 . 

Absorptivity of 
stnicture 


-do- 

-do- 

3. 

Outside Velocity (m/s) 

■ IIIlM 

3.194 

3.194 

•1. 

Inside Velocity (m/s) 


1.2 

1.2 

5. 

No. of Walls 

4 

4 

HHBHHil 

U. 

Volume of Space (m 3 ) 

2638 

1929 

1929 

7. 

No. of Elec.Appl 

10 

2 

1 

8. 

Their Rating (kW) 

20 

10 

■n it 

9. 

No. of Lights 

85 

62 


10. 

Their Rating (kW) 



0.4 

11. 

No. of Occupants 


80 

35 

12. 

East Wall Thick (in) 


0.139 

0.139 

13. 

North Wall Thick (in) 

■KI§gM| 

0.139 

0.139 


11. 

West Wall Thick (m) 

■mgs 

0.139 

0.139 

'15. 

South Wall Thick (m) 

■HKEm 

0.139 

0.139 

1 

10. 

East Wall Area(m' 2 ) 

m 

148 

148 

1 

17. 

North Wall Area(m i ) 

iHKJHM 

47 

47 


18. 

West Wall Area(m 2 ) 

IHfli 

148 

148 

19. 

South Wall Area(m 2 ) 


47 

47 

II 

IH9 

Windows East Wall 

2 

2 

4 

II 

21. 

Windows North Wall 

1 1 

1 

1 

II 

W/msm 

Windows West Wall 

2 

2 

4 

II 

23 . 

Windows South Wall 


1 

1 


mWBM 

Data for Roof 

|| 

in shade 

not shaded 

II 

9S1 

Roof Thickness (in) 


0.174 

0.174 

II 

20. 

Area of Roof(in^) 

im 

533 

533 

J 
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